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Introduction 
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In Newtonian theory of gravitation, there are suprisingly few exact solutions for the proper energy of the homo- 
geneous self-gravitating bodies: a) ellipsoid (including sphere as a particular case) fjj], b) concave bispherical lens 
(including spherical segment as a particular case) (KA), and c) rectangular parallelepiped (including cube as a 
particular case) j|. We present here a method to obtain the formulas for potential energy of homogeneous bodies. 
The method includes a notion of negative density. Here we use this method, first, to check the formulas of KA and, 
second, to get new formula for the homogeneous concavo-convex lens. To this end, we widely used Mathematica M. 



Homogeneous sphere and segment 
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Here we use the method of negative density (MND) to check the formula by KA for potential energy of homogeneous 
spherical segment. 




FIG. 1. Method of negative density. We put into homogeneous sphere, of radius R = OA — OC = OD and density pi, a 
spherical segment of the same radius R, of height h — AB and radius of base a = BC. As a result, we have inhomogeneous 
spherical body comprised of two homogeneous segments with the same radius and base, R and a, respectively. One segment 
has density pi + pi and height h, the other segment has density pi and height 2 R — h = BD. 



In Fig. 1 we show the setting of the problem. First, we have the homogeneous sphere with radius R = OA = OC 
OD and density p\. Gravitational potential, ip, inside the homogeneous sphere, as function of radius r is: 



-2irG Pl (R 2 -r 2 /3); 
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and the gravitational energy of the sphere is: 
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Now we put into the sphere the additional spherical segment with density p2, the same radius R, with height h = AB 
and radius of base a — BC. Note that the total matter density in the region occupied by this segment is p\ + P2- 
Potential energy of interaction of the sphere and the additional segment is obtained by integration of ip * pi over the 
volume of the segment: 
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W int = J <pp 2 dv= J <p[2-Kp 2 r{r-R + h)]dr\ (3) 

V R-h 

we have: 

„. _ Gh 2 ir 2 (h 3 -bh 2 R + 20R 3 ) PlP2 

W int — — . (4J 

According to KA, the gravitational energy of the homogeneous spherical segment with radius R, height h, radius of 
base a and density p 2 is (a = (2Rh — ft, 2 ) 1 / 2 ): 

W segm ent (a, h, R, p 2 ) = -\irGpl^-<K - R) + 8 [3i? 4 a - R 2 a 3 - fa 5 - 6R 4 (R - h) arctan -]]. (5) 

9 2 5 5 a 

If this formula is correct, then the gravitational energy of the homogeneous spherical segment with radius R, height 
2 R— h, base's radius a and density pi is: W a d.se g = W se gment (a, 2 R — h, R, pi). We shall show that this is indeed 
correct. 

We may look at the body in the Fig. 1 in twofold way: as 

a) homogeneous sphere plus homogeneous segment, or as 

b) two homogeneous segments with same a and R but with different heights and densities. 

We can not calculate the total gravitational energy of the body in the case b) as we do not know an external potential 
of spherical segment. However, we can calculate the total gravitational energy of the body in the case a), W a , as a 
sum of three terms: 

the proper gravitational energy of homogeneous sphere, Eq. (2) 

the proper gravitational energy of homogeneous segment, Eq. (5), and 

the potential energy of interaction of sphere and segment, Eq. (4). 

W a = W sphere + W S e 9me „t + W int . (6) 

Now the essence of the negative density method: we take p 2 = —p\\ 

Then, turning to the case b), the body consisted of two segments gives us one (right) "segment" with zero density, and 
the other (left) homogeneous segment with "ordinary" positive density p\. In fact, we have only one (left) segment 
with parameters a, 2 R — h, R, p\. It means that if we take p 2 = —pi, we should get , from Eq. (6), the proper 
potential energy of "left" segment. 

As a result, we have the identity for formula of potential energy of homogeneous spherical segment: 

Gvr 2 (h-R) (-h 2 + 2hR + AR 2 ) 2 p 2 

W se gment{a,2R - h, R, p) - W S egraent{a,h, R, p) = — , (7) 

which is correct if we take into account the evident relation: arctan - + arctan = j . This is a rather rigorous 

check of the formula (||) for potential energy of spherical segment. We conclude that the method of negative density 
is a powerful method of checking the sophisticated formula of gravitational energy of homogeneous spherical segment. 
However our MND may be also applied to check the even more complex formulas by KA for the proper potential energy 
of symmetric and asymmetric concave lenses, and also to obtain the formula for potential energy of concavo-convex 
lens. 
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Homogeneous sphere and asymmetric convex lens 




In Fig. 2 wc show the setting of the problem. First, we have the homogeneous sphere with radius R = OA = 
OC = OD and density p\. Then, we put into this sphere the homogeneous, with density p 2 , asymmetrical bispherical 
lens, comprised of two homogeneous spherical segments: one (right) segment of the same radius R, of height h = AB 
and radius of base a = BC; the other (left) segment of radius Ri > R, of height hi = BE < h = AB and the 
same radius of base a = BC. As a result, we have an inhomogeneous spherical body comprised of two homogeneous 
parts: one asymmetric convex lens (at right) and other asymmetric bispherical concavo-convex lens with surfaces 
having the radii of curvature R (left surface) and —Ri (right surface). The asymmetric convex lens has density 
Pi + pi and central thickness h + hi, the other asymmetric concavo-convex lens has density pi and central thickness 
hi = DE = 1R-h-hi. 

This last figure is of a new kind of homogeneous figures and we are going to calculate the proper gravitational potential 
energy of this lens using MND provided the gravitational energy of asymmetric convex lens is known. 
According to KA (see their Eq. (92)) the gravitational energy of homogeneous asymmetric lens is: 

d = R 2 + Ri - h 2 - hi; a = ^/2R 1 h 1 -h\ = ^2R 2 h 2 - hi; 

Wasl(p 2 , a, Ri,R 2 ) = -^{n[a 2 (h 2 + hi) + R 2 h\ + Rih\\\R\ + R 2 + R 2 h 2 + Rihi - ^-}+ 
■K[hi{Ri - hi)(a 2 + Rihi) - h 2 (R 2 - h 2 )(a 2 + R 2 h 2 )](h 2 -h 1 + R 1 - R 2 )- 

Ti[hl{2R 2 - R 2 h 2 + f ) + h\(2R 2 - Rihi + f )] - MR\h\ + Rfh 2 )- (8) 
^M+M2 [h 2 (3R 2 - h 2 ) + h 2 {ZRi - hi)] - fa 5 - 8a 3 (R 2 + R 2 )- 
^[R 2 (R 2 - h 2 )(R 2 + la 2 ) + R 2 (Ri - hi)(R 2 + la 2 )} + 24a(i?| + Rf)+ 

16i?f[^f - 3(R 2 - /i 2 )]arctan^ + 16 Rf - 3(Ri - hi)] arctan ^}; 

here d is distance between two spheres "generating" asymmetric lens. Potential energy of ASL in the gravitational 
field of homogeneous sphere is (we assume here that p 2 = —pi): 

W intA SL = ^^(Rl +R2- d) 2 

(V + 2 d 3 (Ri + R 2 ) + 5 (Ri - R 2 ) 2 (7R 2 -R 2 )- (9) 
6 d 2 (2 R\ - Ri R 2 + 2 R 2 ) - 2 d (13 R 3 + 9 R 2 R 2 + 9 i?i R 2 - 7 R 3 )^j . 

Not repeating the procedure of previous section we present here only the formula for potential energy of concavo-convex 
lens in the next compact form: 

WcCL = Wsphere + W int ASL ~ W A SL{-pi,CL, Rl,R 2 ); (10) 

see Eqs. (2), (9), (8). 
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Homogeneous sphere and concavo-convex lens 




In this section we calculate the potential energy of asymmetric convex lens provided the potential energy of concavo- 
convex lens is known. The setting of the problem is shown in Fig. 3. Here we attach to homogeneous sphere (at left) 
of radius R± and density p\ , the concavo-convex lens with radii of surface R\ and R 2 and density p 2 . Now we should 
calculate potential energy of CCL in the (external) field of sphere. We have (d is the distance 0\0 2 between centers 
of spheres, R\ — R2 < d < R\-\- R 2 : 

4Gn 2 p lP2 R 1 3 {d-R 1 -2R 2 ) (d-R 1 +R 2 f 

WintCCL = • (11) 

Taking p 2 = pi we get the formula for potential energy of homogeneous figure in Fig. 3 as a sum of three terms (see 
Eqs. (2), (11), (10)): 

W sphere + WintCCL + W C CL, (12) 

which gives the potential energy of ASL coinciding with Wasl(Pi, a, Ri,R 2 ), see Eq. (8). 



Conclusion 



In conclusion, we presented here the method of negative density which is a powerful tool to check the complicated 
formulas by KA for potential energy of homogeneous bispheric lenses. The method allows to obtain new formulas for 
the potential energy of homogeneous bodies as well. Though here we used MND only to study homogeneous bispherical 
lenses, the MND allows to obtain many new formulas, e.g. potential energies of spheres (and ellipsoids) and rectangles 
with spheres, ellipsoids and rectangles laying within (or adjacent to) them. The density p 2 of "additional" bodies may 
be different from density p\ of the "parent" bodies or, according to MND, may be equal to —p\. We leave the study 
of these so called "simple" bodies for next communications. 
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